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ABSTRACT 

Using extremely deep PACS 100- and \60pm Herschel data from the GOODS-Herschel program, we identify 21 infrared bright 
galaxies previously missed in the deepest 24 pm surveys performed by Spitzer [MIPS. These MIPS dropouts are predominantly found 
in two redshift bins, centred at z ~0.4 and ~1.3. Their S 100/^24 flux density ratios are similar to those of local (ultra-) luminous 
infrared galaxies (LIRGs and ULIRGs), whose silicate absorption features at 18//m (at z ~ 0.4) and 9.7 pm (at z ~ 1.3) are shifted 
into the 24 pm MIPS band at these redshifts. The high-z sub-sample consists of 1 1 infrared luminous sources, accounting for ~ 2% 
of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0 < z < 1.7 that are 
missed by the 24 pm surveys. We find that a S100/S24 > 43 colour cut selects galaxies with a redshift distribution similar to that of 
the MIPS dropouts and when combined with a second colour cut, Si 6 /S% > 4, isolates sources at 1.0 < z < 1.7. We show that these 
sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be 
compact starbursts with moderate/strong 9.7 pm silicate absorption features in their mid-IR spectra. Herschel data reveal that their 
infrared luminosities extrapolated from the 24 pm flux density are underestimated, on average, by a factor of ~ 3. These silicate break 
galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density 
of 2.0 x 10~ 5 Mpc -3 . Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 
250- and 350 pm sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the 
major 24 pm extragalactic surveys have been overestimated. 
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1. Introduction 

Accurate measurement of star formation rates (SFR) is a key 
ingredient for studying galaxy evolution and deriving the census 
of the star formation activity, both in the distant and in the local 
universe. To this end, it has been shown that the contribution of 
luminous infrared galaxies (with infrared luminosities Lir > 10 11 
L©) to the star formation density is progressively rising as we 
look back in cosmic time, at least up to z ~ 2. Indeed, although 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



they were found to be rare in the local Universe and to account 
for only ~ 5% of the total infrared energy emitted by galaxies 
at z ~ (Soifer et al. 1991, Kim & Sanders 1998), LIRGs and 
ULIRGs (Lir > 10 12 L Q ), dominate the SFR density at z ~ 1-2, 
accounting for 70% of the star formation activity at these epochs 
(Papovich et al. 2004, Le Floc'h et al. 2005, Caputi et al. 2007). 

The study of infrared sources was greatly facilitated by the 
advent of the Spitzer Space Telescope (Spitzer Werner et al. 
2004). Extragalactic surveys carried out using the MIPS 24 pm 
band on-board Spitzer confirmed the strong evolution of these 
sources first indicated by infrared and submillimeter observa- 
tions using ISO and SCUBA, respectively (Blain et al. 1999a, 
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Elbaz et al. 1999, Serjeant et al. 2001, Dole et al. 2001). Such 
surveys are believed to detect the bulk of the dusty star forming 
galaxies up to z ~2. However there are two important caveats. 
The first is that the conversion of 24 jim flux densities to total 
Lir (and therefore SFR) is subject to large uncertainties, as it re- 
lies on extrapolations that strongly depend on the assumed spec- 
tral energy distribution (SED) libraries (Chary & Elbaz 2001; 
Lagache et al. 2003; Dale & Helou 2002). The second comes 
from the prominent emission and absorption features between 
3- and 19 yum in the spectra of star forming galaxies that parade 
through the 24 yum band at various redshifts. 

A large number of studies using the Infrared Spectrograph 
(IRS, Houck et al. 2004) have revealed that the vast majority 
of local LIRGs and ULIRGs exhibit a broad silicate absorption 
feature centred at 9.7 yum with silicate optical depths ranging 
from T9.7 ~ 0.4 to r 9 . 7 > 4.2 (Brandl et al. 2006; Armus et al. 
2007; Pereira-Santaella et al. 2010). Furthermore, using a sam- 
ple of local ULIRGs, Desai et al. (2007) found strong PAHs 
and prominent silicate absorption in the H II and LINER sources 
and weak PAHs and silicate absorption in Seyferts, suggesting 
that ULIRGs with strong PAHs but weak silicate absorption are 
rare. Similar features have been observed in the mid-IR spec- 
tra of high-z galaxies (Higdon et al. 2004, Houck et al. 2005). 
For example, Menendez-Delmestre et al. (2009) and Farrah et al. 
(2008) report a median T9.7 ~ 0.31 for a sample of submillimetre 
galaxies (SMGs) and IRAC selected ULIRGs respectively, while 
Sajina et al. (2007) found deeper silicate absorption features (T9.7 
> 1.1) in a sample of z ~ 2 radio-loud galaxies. Although for 
LIRGs and ULIRGs the identification and measurement of the 
silicate optical depths is straightforward, this is not the case for 
normal galaxies (Li R < 10 10 L ), as for the latter, it is difficult 
to discriminate between moderate PAH emission superimposed 
on a silicate-absorbed continuum and strong PAH features with 
a relatively weak underlying continuum (Smith et al. 2007). 

Whatever its origin, the existence of this broad dip in the 
mid-IR spectra of star forming galaxies 10 jim could be of partic- 
ular importance for galaxies in the redshift range of 1 < z < 1.8. 
At these redshifts the 24 yum filter samples this part of the spec- 
trum and sources with such features would appear faint at 24 jim 
or even be undetected in this band (depending on the depth of the 
24//m data). A second broad dip that is common in the spectra 
of star forming galaxies is caused by another silicate absorption 
feature at 1 8 yum. This feature would have a similar effect for 
sources at 0.2 < z < 0.6. 

The impact of these features on the mid-IR colours as a func- 
tion of redshift were presented in detail by Takagi & Pearson 
(2005), who predicted a population of infrared luminous galax- 
ies at z ~ 1.5 which, due to strong absorption at 9.7 yum, are not 
detected in the 24/im band. Subsequent studies that focussed 
on the search for such silicate absorbed systems employed the 
16//m IRS peak-up image. Kasliwal et al. (2005) suggested that 
such objects account for more than half of all the sources at z ~ 
1-2 predicted by various models. It has also been proposed that 
the mid-IR colour anomalies introduced by the silicate absorp- 
tion feature can serve as a redshift indicator for dusty infrared 
luminous galaxies at z -1.5 (Charmandaris et al. 2004, Teplitz et 
al. 201 1, Armus et al. 2007). Similar claims have also been pre- 
sented by Pearson et al. (2010), using the AKARI IRC L18W 
to MIPS24 band colour. These studies raised concerns about a 
possible bias introduced by the 24 jim selection, in the sense that 
a significant fraction of z < 2 LIRGs and ULIRGs could remain 
undetected in 24 /mi surveys. However, with little or no infor- 
mation about the far-IR part of the spectrum, these studies were 



subject to large extrapolation and hence suffered from large un- 
certainties. 

With the successful launch of the Herschel Space 
Observatory (Herschel, Pilbratt et al. 2010), we now have ac- 
cess to wavelengths that directly probe the peak of the far- 
IR emission of high-z galaxies and are in a position to mea- 
sure with unprecedented accuracy their bolometric output. Deep 
Herschel extragalactic surveys can be used to determine the ac- 
curacy of our extrapolations of the far-IR properties of high-z 
galaxies as well as test previous claims that 24 jim surveys miss 
a population of z < 2 LIRGs and ULIRGs (24 yum dropouts). In 
this paper, we use the deepest Herschel observations to date, 
as part of the GOODS-Herschel (GOODS-H) program (PI D. 
Elbaz), covering both the north and the south part of the GOODS 
fields (GOODS-N and GOODS-S respectively), (Dickinson et 
al. 2003, Giavalisco et al. 2004), to search for such sources. In 
Section 2 we present the Herschel data, introduce the GOODS- 
H sample of galaxies, and identify 24/jm dropout sources, i.e., 
sources detected in the PACS bands but not at 24 jim. In Section 
3 we investigate the properties of this population, while in 
Section 4 we extend our study to the whole GOODS-H sam- 
ple. Finally in Section 5 we provide estimates of the fraction of 
z <2 MIPS dropout sources as a function of the 24yum, lOO/mi 
and 160yum sensitivity limits and summarize our results. 

2. Herschel data and sample selection 

Herschel observations were obtained as part of the open time key 
program GOODS-H (PI D.Elbaz). The full 10' x 16' GOODS- 
N field was imaged with the PACS (Poglitsch et al. 2010) and 
SPIRE (Griffin et al. 2010) instruments at 100, 160 yum (PACS) 
and 250, 350, 500 /mi (SPIRE). The total observing time was 
124.6 hours (~ 2.5h / sky position) and 31.1 hours for PACS 
and SPIRE respectively. Similarly a7x7 part of the GOODS- 
S field was observed by PACS over a total of 264 hours (~ 
15h/ sky position). Observations of both fields were carried out 
by adopting the intermediate speed (20" s -1 ) scan-map mode. 
Both PACS and SPIRE data were processed through the standard 
Herschel reduction pipeline, version 6.0.3, within the HCSS en- 
vironment. Additionally, we employed custom procedures aimed 
at removing of interference patterns, tracking anomalies, re- 
centering positional offsets, and mapping. A full description of 
the data reduction procedures will be given in a companion paper 
(Leiton et al. 201 1 in prep). 

2.1. Prior based Source extraction; The GOODS-H sample 

Given the large beam size of the Herschel bands, (FWHM ~ 
6.7", 11.2",18.0", 25.0", 36.0" for PACS 100- and 160 ^m and 
SPIRE 250, 350 and 500 yum), a common approach to perform- 
ing source extraction has been a guided extraction using pri- 
ors. Here we will give a brief summary of the procedure as 
an extensive description of the method is given in Elbaz et al. 
(2011). Source extraction and photometry were obtained from 
point source fitting at prior positions defined by 24 jim sources 
with fluxes brighter than S24 ~ 20 yuJy for the 100 jim maps and 
down to S 24 ~ 30 yuJy for the 160yum and 250 /mi maps. For 
the other two SPIRE bands, a secondary criterion was needed, 
as the 24 jim sources were far too numerous and would lead to 
an over-deblending of the actual sources. Hence, only sources 
with S/N > 2 at 250 jim were considered as priors for the longer 
wavelength SPIRE bands. This choice was optimized by Monte 
Carlo (MC) simulations to avoid artificial over-deblending of a 
source, but also to give clean residual maps. Flux uncertainties 
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Fig.l. ACS V-band (5" x 5"), IRAC 3.6/mi (20" x 20"), MIPS 24 //m (30" x 30"), and PACS 100-, 160/rni (30" x 30") cut-out images of the 
24 //m dropout sources from our sample. The red circles are centred at the IRAC 3.6 //m positions of the sources and their diameter corresponds to 
the FWHM at each band. The size of each image is denoted on the top of each column. 



were based on local estimates of the background noise at the 
position of the sources from residual images produced after sub- 
tracting detected sources, while global noise estimates for the 
maps were derived from Monte Carlo simulations. This used ar- 
tificial sources injected into the Herschel maps and source ex- 
traction performed in the same manner as for the real sources. 
The dispersion between input and recovered fluxes provides a se- 
cure estimate of the completeness and the noise properties of the 
map. The two noise estimates were found to be in good agree- 
ment. 

To construct the GOODS-H sample, we considered sources 
with flux densities down to 3cr in the PACS bands, i.e. 1.0 and 
2.6 mJy (0.7 and 2.6 mJy) at 100 and 160 /mi in GOODS-N 



(GOODS-S). For the GOODS-N sample, where SPIRE data are 
available we also considered sources down to the 5<x detection 
limit, i.e. 6.3, 7.1, and 15.0 mJy at 250, 350 and 500 /mi, respec- 
tively. The choice of a higher S/N cut for the SPIRE catalogues 
was dictated by the larger beam size and the confusion noise 
that significantly affects the SPIRE observations (for more de- 
tails see Elbaz et al. 201 1). Herschel catalogs were then matched 
with the existing multi- wavelength data of the GOODS team to 
create a multi-band merged catalogue of GOODS galaxies in- 
cluding HSTACS BViz, J,K, IRAC, Spitzer MIPS 24 and 70 /mi, 
Herschel PACS 100 and 160 /mi, and Herschel SPIRE 250, 350 
and 500 /mi. Among our sources, 65% have secure spectroscopic 
redshifts, while for the rest, we use the reliable compilation of 
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photometric redshifts by Le Borgne et al (2009). Hereafter we 
will refer to this sample as the GOODS-H sample. 



2.2. Blind source extraction; The 24 jim dropout sample 

Since the main aim of this work is to investigate whether 
Herschel observations reveal a population of galaxies that were 
previously missed by 24 /mi surveys, we also performed blind 
source extraction in the two PACS bands using Starfinder, a 
point spread function (PSF) fitting code (Diolaiti et al. 2000). 
We first extracted PSF profiles from the final science maps that 
were used to perform source extraction. Aperture corrections 
were derived based on calibration observation of the asteroid 
Vesta, while the flux uncertainties were derived based both on 
the error maps and Monte Carlo simulations, as described above. 
Monte Carlo simulations were also employed to obtain the level 
of completeness and the fraction of spurious sources. Both the 
derived fluxes and the noise properties of the maps are in good 
agreement with those obtained by the prior based source extrac- 
tion. Finally, a critical parameter in Starfinder is the correlation 
threshold (ct), a measure of the similarity between the PSF used 
for source extraction and the profile of the extracted source, with 
ct=l corresponding to identical profiles. MC simulations indi- 
cate that high ct values result in catalogs immune to spurious 
detections but with lower completeness levels. Similarly lower 
ct values correspond to higher completeness but also to higher 
fractions of spurious sources. For our blind catalogs we consider 
sources with flux densities above the 4<x detection limit at each 
band and ct > 0.67, for which the fraction of spurious sources is 
< 3%. 

In principle, our aim was to find sources detected in either 
of the PACS bands but undetected at 24 /mi. Therefore, we first 
matched the PACS 100- and 160yum blind catalogs with the 
MIPS24yum sample down to S24 ~20 /Jy (3cr), i.e. the one that 
served as a pool for the prior based source extraction, starting 
from the longest wavelength available and using search radii of 
7" and 11" respectively. Sources with 24yum counterparts were 
omitted while the rest were matched to the IRAC 3.6 yum cata- 
logue and subsequently to the master GOODS multi-bands cata- 
logue described in the previous section. We also performed pho- 
tometry in the 24 yum maps at the position of the PACS sources 
to ensure that there is no 24 yum source at this position, possibly 
missing from the 24 yum catalogue. We calculated the corrected 
Poissonian probability that an association of a PACS source 
within the search radius is a chance coincidence (see Downes et 
al. 1986) and all sources were found to have a robust (p < 0.05) 
3.6 fim counterpart. We also note that flux boosting due to in- 
sufficient de-blending, which is the main caveat of blind source 
extraction, should not, by definition, be an issue for the 24 jim 
dropout sample. All sources were also inspected by eye and a 
quality flag was attributed to them. In particular, sources with 
multiple IRAC counterparts within the PACS beam and sources 
close to bright objects in the PACS bands were flagged as low 
quality sources. The final sample consists of 21 MIPS dropout 
sources, all detected at lOOyurn and two at 160yum, account- 
ing for ~ 2% of the total sources detected in the PACS bands. 
Hence, we find that even at the confusion limit of the 100- and 
160 yum passbands (0.7- and 2.6 mJy at a 3<x level), 98 % of 
the Herschel sources possess a robust 24 jim counterpart brighter 
than 20 yuJy. The 21 MIPS dropout sources are shown in Fig. 1 
where we present cut-out images at several bands. 



3. 24yum dropout sources 

The small number of 24 yum dropouts indicates that the vast ma- 
jority of PACS sources do have a 24 yum counterpart. In other 
words, the "normal" SED behaviour of galaxies in the GOODS 
sample is the one where the relative sensitivity at 24 yum over- 
powers that of the PACS bands. In Fig. 2 (left) we plot the flux 
density at 24yum over that at lOOyurn for the whole GOODS-H 
sample, as well as for the dropout sources and see that the latter 
depart from the general trend and are relatively faint at 100 /mi. 
We also note that similarly to the MIPS dropouts, some sources 
with 24yum detection tend to exhibit redder S 100/S 24 colours 
than the bulk of the population while they span a wide range of 
S24. Here we will study the origin of the departure of the MIPS 
dropouts from the bulk of the GOODS-H sample and our investi- 
gation will be driven by their property that intrigued our interest 
in the first place i.e. their unusual S 100/^24 colour. 

3.1. Far-IR properties 

The total infrared luminosity (Lir = Lg-iooOyum) °f galaxies in the 
sample was determined from the lOOyurn flux density using the 
templates of (Chary & Elbaz 2001, CE01) and Dale & Helou 
(2002) (DH02). Despite the lack of data points at longer wave- 
lengths, we note that the monochromatic derivations of total IR 
luminosities from the lOOyurn flux density tend to be robust up 
to z ~ 1.5 (Elbaz et al. 2010). For the two sources with 160 yum 
detection Lir was determined from the best fit of the two PACS 
points, using the whole library of SED templates from CE01 in- 
dependently of their luminosity (i.e. allowing normalization of 
all SEDs to the observations) as well as the DH02 templates. 
The derived luminosities range from 8 x 10 9 to 2 x 10 12 L©, with 
11 of the MIPS dropout sources having L IR >10 n L©, belong- 
ing to the class of luminous infrared galaxies. To illustrate the 
completeness of our sample in terms of Lir, we plot the Lir for 
the GOODS-H sample, as well as for the drop-out sources as 
a function of redshift, along with the corresponding detection 
limits at 100 yum and 24yum (Fig. 2 middle). Furthermore, In 
Fig. 2 (right), we also plot the S 100/S 24 colour as a function 
of Lir as derived from the Herschel data, both for the whole 
GOODS sample as well as for the MIPS dropouts (crosses with 
arrows). For the latter we compute lower limits to S ioo/S 24 as- 
suming the 3cr detection limit of the 24 jim maps. The points are 
also colour-coded based on their redshift. We see that 24 fim 
dropouts have Lir values similar to that of the whole GOODS- 
H sample (for a given redshift), while they exhibit significantly 
higher S 100/^24 colours for the whole range of luminosities. As 
we discussed above, Lir scales with S 100 . Hence, given the rich- 
ness of features in the rest frame MIR emission (i.e. Armus et 
al. 2007), it is more likely that the dropouts have a suppressed 
S 24 emission rather than an excess at S 100 when compared to the 
rest of the GOODS-H sample. In what follows we investigate the 
origin of this S 24 deficit. 

3.2. S 100/S24 colour and redshift distribution 

We wish to examine whether the MIPS dropout sources tend to 
be found at specific redshifts. We note that nine sources have 
spectroscopic redshift, while for the rest we adopt the photomet- 
ric redshifts derived by the GOODS team. It appears that 24 jim 
dropouts are distributed in two redshift bins, one centred at z ~ 
0.4 and one at z ~ 1.3 (Fig. 3). This bimodality is verified by 
a KMM test (Ashman et al. 1994) at a 4.3 cr confidence level. 
Therefore, it seems that our sample is populated by low-z (0.2 < 
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Fig. 2. Left: S 24 vs S 100 flux densities for the whole GOODS-H sample (black circles) as well as for the 24/im dropout sources (red arrows). 
For the dropouts we consider a 3 cr upper limit for the S 24. Dropouts as well as some 24/um detected sources tend to depart from the bulk of the 
GOODS-H population, exhibiting redder S 100/^24 colours. The cyan line corresponds to S 100/^24 = 43. As discussed latter in the paper, sources 
with S iqq/S 24 > 43 are classified as silicate-break galaxies. Middle: Detection limits as a function of redshift for the GOODS-N and GOODS-S 
PACS 100//m and MIPS 24yt/m observations. Red squares correspond to the drop-out sources. Right: S100/S24 as a function of L IR as derived 
by Herschel for the whole GOODS-H sample (circles) and lower limits for the 24yi/m dropouts (arrows). Both samples are colour coded based 
on their redshift. Sources with a black cross are AGNs based on their X-ray emission. Filled symbols denote sources in the GOODS-H sample 
with spectroscopic redshifts while open symbols sources with photometric redshift. Similarly, yellow circles on top of the arrows indicate that a 
spectroscopic redshift is available for that 24 fim dropout source. The horizontal black dashed line corresponds to S 100/^24 = 43. 
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Fig. 3. Redshift distribution of sources in the MIPS dropout sample. A 
KMM test suggests a bimodal distribution centred at z ~ 0.4 and z ~ 
1.3. Blue shadowed area corresponds to the distribution of sources with 
spectroscopic redshift. 



z < 0.6) and high-z (0.9 < z < 1.7) sources. Similarly to Fig. 2, 
we now plot the S 100/ S 24 colour as a function of redshift (Fig. 4). 
In this figure we also overplot the S 100/^24 colour as a function 
of redshift for a number of local LIRGs/ULIRGs based on their 
observed SED as constructed by IRS observations of their mid- 
IR spectrum (Armus et al. 2007) and IRAS observations of their 
far-IR emission (Rieke et al. 2009). The observed templates were 
chosen to span a wide range of Lir and silicate optical depths 
(T9.7) and their full SEDs are presented in Fig. 5, in order of in- 
creasing T9.7 . The T9.7 measurements are adopted from Armus et 
al. (2007) (Arp200, NGC 22491), Pereira-Santaella et al. (2010) 
(ESO 320-G030, NGC 2369, Zw 049.057, NGC 3256) and da 
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Fig. 4. S 100/S24 as a function of redshift for the whole GOODS-H sam- 
ple (circles) and lower limits for the 24 11m dropouts (arrows). Both 
samples are colour coded based on their L IR . Sources with a black cross 
are AGNs based on their X-ray emission. Filled symbols denote sources 
in the GOODS-H sample with spectroscopic redshift, while open sym- 
bols are sources with photometric redshift. Similarly, black circles on 
top of the arrows indicate that a spectroscopic redshift is available for 
that MIPS dropout source. Solid lines correspond to different observed 
SEDs of local LIRGs/ULIRGs (see Fig. 5) and horizontal black dashed 
line to S 100/^24 = 43. The pop panel shows the SED of Arp220 at vari- 
ous redshifts along with the MIPS 24- and PACS 100 jjim bands. 
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Cunha et al. (2010) (IRAS 17128). On the top of Fig. 4 we also 
show the SED of Arp220, along with the 24- and 100 /mi bands 
at several redshifts. It is evident that the S 100/^24 colour of the 
templates varies significantly as a function of redshift, mainly 
due to the presence of the silicate absorption features at 9.7- 
and 18 /mi that enter the MIPS 24 /mi filter at z ~ 1.4 and ~ 
0.3 respectively. We also note a wide range of S 100/S 24 colours 
for a given redshift, indicative of different amounts of extinction 
as well as different dust temperatures. Indeed, Fig. 5 suggests a 
weak trend for SEDs of sources with deeper silicate absorption 
to peak at shorter wavelengths. Therefore, the S 100/^24 values 
increase because S 24 is suppressed by the silicate absorption fea- 
tures, but are also further elevated due to higher dust (big grain) 
temperatures and hence S 100 values. 

Looking back at Fig. 4, we see that 24 /mi dropouts, although 
they have higher S ioo/S 24 values when compared to the whole 
sample, exhibit colours that are consistent with those of local 
star formation dominated ULIRGs/LIRGs. The fact that the red- 
shift distribution of 24 /mi dropouts peaks at redshifts where the 
9.7- and 18 /mi silicate absorption features enter the 24 /mi band 
indicates that these sources could have moderate/strong silicate 
absorption features. On the other hand, the high S 100/S 24 ratio 
could also result from lower levels of observed dust continuum 
emission at 24 /mi, but as we will see later that our data dis- 
favour this scenario. In any case, all 24 /mi dropouts fall within 
the envelope defined by the templates, so they form an extreme 
rather than an extraordinary population of star forming galaxies. 
We also note that none of the drop-out sources have an X-ray 
detection or meet the criteria for a power-law AGN. 

Previous studies have attempted to identify infrared lumi- 
nous sources that are undetected at 24 /mi by combining 16 /mi 
IRS peak-up imaging with 24 /mi data (Kasliwal et al. 2005, 
Teplitz et al. 201 1) and following an approach similar to the one 
presented in this work. They searched for sources at z ~ 1.3, that 
are faint at 24 /mi, due to the shift of the silicate absorption fea- 
tures into the MIPS band, but bright at 16 /mi due to 7.7 /mi PAH 
emission, concluding that sources with S 16/S24 > 1.2 tend to be 
found at z >1.1. They also investigated the Sie/S 24 colours of 
several local LIRGs/ULIRGs and report an average ratio of ~1- 
2 for galaxies with strong silicate absorption features. From our 
24 fim dropout sample, none of the sources are detected at 16 /mi 
down to the 3cr detection limit (S i6 ~ 40 /Jy for GOODS-N and 
~ 65 /iJy for GOODS-S, Teplitz et al. 201 1). Apart from the fact 
that some of our sources are outside the area covered by the IRS 
peak-up image in the GOODS fields, the non-detection of the 
rest of the sources at 16 /mi is somewhat expected from the dis- 
cussion above. Even if we adopt a S24 = 20 //Jy for all sources 
in the sample and assume a S 16/S24 = 2.0 (Arp220 case, Armus 
et al. 2007) then our sources should only be marginally detected 
at the depth of the 16 /mi GOODS maps. 

3.3. Az~ 1.68 ULIRG, missed by MIPS 

Although the study of the properties of individual 24 /mi dropout 
sources is beyond the scope of this study, here we wish to have a 
closer look at a specific source for which, Herscheldata. indicate 
an infrared luminosity Lir > 10 12 L . This is the only source 
in our high-z sample with detection at 100 /mi (Sioo = 5.2 ± 
0.42 mJy), 160 /mi (S m = 6.9 ± l.lmJy) and at 250/mi. The 
source also has also a ~ 6<x detection at the GOODS-N VLA 
1.4GHz map (5 1.4 = 32.1 ±5.1 //Jy). Unfortunately, the source 
is out of the area covered by IRS 16 /mi peak-up imaging. Cut 
out images of this galaxy are shown in the sixth row of Fig. 1. 
To derive the far-IR properties of the source we fit the observed 
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Herschel points with CE01 models, but also with a range of ob- 
served SEDs of local ULIRGs, described in the previous section. 
The observed multi-band photometry, together with the best fit 
SEDs are shown in Fig. 6. 

The derived infrared luminosity of the source is Lir = 
1.3-1.6 x 10 12 L© (depending on the assumed SED), but it is 
evident that only templates with a strong silicate absorption fea- 
ture at 9.7 /mi can reproduce the non-detection down to 20 //Jy 
at 24 /mi. Fitting the optical part of the SED with BC03 models, 
yields a stellar mass of 8.9 (±0.5) x 10 10 M Q , while based on 
the J3 slope of the UV spectrum we derive a reddening E(B-V) 
= 2.1, consistent with a heavily obscured source. Assuming a 
Salpeter IMF and the Kennicutt (1998) relation we convert the 
Herschel-based Lir to star formation rate and derive SFR ~ 280 
M© yr" 1 . Similarly, we convert the radio flux to Lir (Condon 
1992) and subsequently to SFR, finding that the two estimates 
are in perfect agreement (S FR ra & - 300 ± 32 M© yr" 1 ). The 
source is not detected in X-rays and exhibits the 1.6 /mi stel- 
lar bump, indicating that there is no near-to-mid-IR evidence 
for the presence of an AGN. The above analysis suggests that 
this is a heavily obscured starburst galaxy with no signs of AGN 
activity. Although unique in our sample, this galaxy raises the 
question of how many such objects we might have missed in 
the pre-Herschel era and what would be their contribution to the 
cosmic SFR density. A detailed discussion on this point will be 
presented in section 5. 



4. Silicate-break galaxies 

So far, we have demonstrated that Herschel data have revealed a 
small but interesting class of infrared luminous galaxies, which 
are undetected at 24 /mi. The main characteristic of these sources 
is their atypically red S 100/^24 colour. Furthermore, Fig. 2 (left), 
reveals that sources with S ioo/S 24 colours similar to those of the 
dropouts can be found among sources with 24 /mi detection. In 
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Fig. 6. The SED of a dropout at z = 1.68, detected in both the 100- and 
160yum bands, although missed at 24 iim. The optical part is overlaid 
with the best fit BC03 model (green line) and the infrared part with a 
range of observed and model SEDs. The red square denotes the upper 
limit at lAfim. Cut-out images of the source are shown in Fig. 1, 6th 
row. 
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Fig. 7. Top: Redshift distribution of sources in GOODS-N with PACS 
100/zm and MIPS 24/im detections. Bottom: Redshift distribution of 
sources in GOODS-N with S 100/^24 > 43. Sources are clustered around 
z = 0.4 and z =1 .3. The blue histogram corresponds to the redshift 
distribution of MIPS dropouts. 



what follows we extend our study to the whole GOODS-H sam- 
ple, searching for such objects. 

To select them, we adopt a cut off, S 100/^24 > 43, that corre- 
sponds to the bluest lower limit among the MIPS dropout sample 
and the colour of the local LIRG ESO320 shifted at z = 1.4. In 
Fig. 7 we plot the redshift distribution of sources with S 100/S 24 > 
43 (excluding the dropouts), finding that it bears a remarkable 
resemblance to that of the 24 /mi dropouts in Fig. 3. Indeed, a 
K-S test reveals that there is no significant difference between 
the two samples with a p value of 0.61. In contrast, we find that 
the sample has a redshift distribution quite different than that of 
the whole GOODS-N sample (Fig. 7, top panel) at a confidence 
level of >97.5%. Similar to the dropout sample, sources with 
S 100/S24 > 43 exhibit a bimodal redshift distribution at the 3. So- 
level, centred at z ~ 0.48 and ~ 1.3. Therefore, in what follows, 
we suggest that the S 100/ S 24 colour could serve as redshift indi- 
cator. 

4.1. S 100/S24 colour as a redshift indicator. 

In the pve-Herschel era, there were several attempts to use 
anomalous MIR colours as a crude redshift indicator. These stud- 
ies mainly employed the 16yum IRS peak-up imaging and pro- 
posed that a blue S 16/S24 colour would peak in a narrow redshift 
bin, 1.0 < z < 1.8 (Armus et al. 2007, Teplitz et al. 201 1), as well 
as being useful for selecting infrared luminous galaxies in this 
redshift bin that were undetected by the 24 yum band. The main 
idea behind this criterion is that at these redshifts the 9.7 yum sil- 
icate absorption feature enters in the 24 /mi passband, while the 
strong PAH emission features at 6.2- and 7.7 yum are shifted into 
the 16//m band, producing a distinctive blue S 16/S24 colour. For 
example, Kasliwal et al. (2005) used a Sie/S 24 > 1.2 ratio to 
select objects at 1.1 < z < 1.6, and called them "silicate-break 



galaxies", attributing their blue S 16/S24 colour to the existence of 
a strong silicate absorption feature at 9.7 yum. This approach had 
two main caveats. First, the blue S 16/S24 colour is not necessar- 
ily produced by a silicate absorption feature, as similar colours 
can appear for sources with strong PAHs and low dust contin- 
uum. Second, many objects at lower redshifts fall within the 
same colour cut (Teplitz et al. 2005). A more recent attempt by 
Teplitz et al. (2011), reports that a higher ratio (S i6 /S 2 4 > 1.4), 
would eliminate many but not all of the low-z interlopers. 

Here, we face the same situation. The cut in the S 100/^24 ra- 
tio that we have adopted selects sources in two redshifts bins. 
In order to reject the low-z sources in our sample, simply in- 
creasing the S 100/ S 24 cut is not useful as apart from the low-z, 
we also miss many high-z sources and the selection is still not 
pure enough. Alternatively, we can employ a second colour cri- 
terion, based on the 16/mi and 8 jim flux densities. In Fig. 8 we 
plot the S 100/^24 vs S 16/58 colour - colour diagram for sources 
that have a 16- and 8 /im detection. We see that if a ratio cut of 
S 16/ S 8 > 4 is combined with a S 100/^24 cut > 43, then we suc- 
cessfully reject all low-z interlopers, while selecting sources in 
the 1 .0 < z < 1 .7 redshift bin. The selection of S \e/S g as a second 
colour criterion was driven by the fact that at 1.0 < z < 1.7 the 
16//m band probes the PAH complex (6.2 - 7.7/mi), boosting 
the value of S 16 while for the low redshift sources only traces 
emission from a warm dust continuum at > 10 yum. We note that 
none of our sources is classified as an AGN, based on their X- 
ray emission, their optical spectra or their mid-IR colours (i.e. 
power law AGNs). It therefore seems that we have found a way 
to select star-forming high-z galaxies in a narrow redshift bin. In 
what follows we argue that these sources are compact starbursts 
with moderate/strong silicate features in their MIR spectrum. 
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Fig. 9. Distribution of the S 100/^24 colour among sources in the 
GOODS-H sample with spectroscopic redshift at 1.0 < z < 1.7. Green 
line and green shadowed area indicates the Gaussian fit to the distri- 
bution and the area beneath it. Coloured bars indicate the range of the 
S 100/^24 colour of several observed SEDs of local LIRGs and ULIRGs 
at this redshift bin. The silicate optical depths of the local sources, as de- 
rived by IRS spectroscopy are also overlaid. The black arrow indicates 
the position of the high-z 24 jim dropout sample. 



4.2. Evidence for silicate absorption 

We have already discussed in the Introduction, that for normal 
galaxies, i.e. those with L IR <10 10 L , discriminating between 
moderate PAH emission superposed on a silicate-absorbed con- 
tinuum and strong PAH features with a relatively weak under- 
lying continuum is a difficult task (e.g. Smith et al. 2007), even 
when high quality mid-IR spectra are available. A typical ex- 
ample is M82, for which Sturm et al. (2000) suggests that there 
is no silicate absorption since the lOyum dip can be reproduced 
by a superposition of strong PAHs and VSG continuum. On the 
other hand, things are more straightforward for more luminous 
infrared sources, where there is clear evidence for the existence 
of a wide range of silicate optical depths, both for local and high- 
z LIRGs and ULIRGs (Armus et al. 2007, Pereira-Santaella et al. 
2010). Unlike the sample based on the blue S\e/S 24 colour, for 
which the Lir estimates were based on large/uncertain extrap- 
olations, our study benefits from more robust Herschel-based 
Lir estimates. The fact that all sources with 5 100/^24 > 43 at 1.0 
< z < 1.7 in our sample (including the high-z 24yum dropouts) 
have L IR > 10 11 L , coupled with the lack of observational ev- 
idence of sources with similar luminosities and low dust con- 
tinuum emission, is a first hint that these sources have a silicate 
absorption feature at 9.7 /mi. 

We have already demonstrated in Fig. 4, that these sources 
have 5 100/ 5 24 colours consistent with those of local templates 
of LIRGs and ULIRGs with moderate/high silicate strength. 
To further investigate this, we plot the distribution of the 
S100/S24 colour of the GOODS-H (GOODS-N and GOODS- 
S) sources with spectroscopic redshift at 1.0 < z < 1.7, along 
with the range of the 5 100/ 5 24 colour at this redshift range for 
5 LIRG/ULIRG observed templates with a wide range of T9.7 
values (Fig. 9). In practice, this plot examines the detailed distri- 
bution of galaxies in the 1 .0 < z < 1 .7 region of Fig. 4. Fitting the 



colour distribution with a Gaussian indicates a clear excess in the 
red tail due to sources with 5 100/^24 > 43. It is crucial to stress 
that this excess is absent in other redshifts bins. Furthermore, we 
see that based on the template SEDs, sources with higher opti- 
cal depth exhibit redder 5 100/^24 colours in this redshift bin. In 
other words, it seems that the 5 ioo/5 24 colour could serve as a 
rough indicator of the silicate strength. Our sources, i.e. those 
with 5 ioo/5 24 > 43, have colours consistent with 0.7 < T9.7 < 
1.9, while it appears that none of them have a silicate feature as 
strong as that of Arp220. We note however, that we cannot rule 
out a source in our sample with stronger silicate but colder 7d us t- 
An alternative explanation of the high 5ioo/524 ratios 
of these sources would be enhanced 100//m emission, sim- 
ply because the sources are warmer, and hence their SEDs 
peak at shorter wavelengths. To test this, we investigate the 
5 ioo/5 24 colour as a function of dust temperature, as indicated 
by the 5 i6o/5 100 ratio. In Fig. 10, we present this colour-colour 
diagram for sources at 1.0 < z < 1.7, which are detected at 160- 
and 100/zm. It appears that the adopted 5ioo/524 cut does not 
introduce a bias towards warmer sources (i.e. sources with low 
5 i6o/5 100 colours), as it selects sources with a wide range of dust 
temperatures similar to that found for the whole GOODS-H sam- 
ple in this redshift bin, suggesting that an enhanced 5 100 is not 
the main reason for their high 5ioo/524 ratios. We note that a 
large dispersion of dust temperature of high-z galaxies has re- 
cently been demonstrated by Hwang et al. (2010) and Magdis et 
al. (2010d). 

One of the sources that meets the selection criteria is the 
well studied sub-millimetre source at z = 1.21, GN26 for which 
Pope et al. (2008) have published an IRS spectrum. The detec- 
tion of this source in the PACS bands has also been discussed 
by Dannerbauer et al. (2010) and Magnelli et al. (2010) while 
Frayer et al. (2008) has reported the detection of CO(2— > 1) emis- 
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Fig. 10. *Si6o/^ioo vs S100/S24 colour - colour diagram for sources in 
the GOODS-H sample at 1.0 < z < 1.7. Sources with S 100/^24 > 43 and 
S m/Ss > 4 are filled with orange colour. This plot demonstrates that our 
selection is not biased towards warmer sources. 



sion. In Fig. 11, we show the full SED of the source, overlaid 
with the best CE01 template for the far-IR part and the observed 
IRS spectrum for the mid-IR part of the SED. The IRS spec- 
trum of GN26 is also presented separately in an inset panel. This 
source has S 100/^24 ~ 70, and according to Fig. 9, it should have 
a silicate optical depth of ~ 1 . Although the S/N of the IRS spec- 
trum does not allow for a robust measurement of T9.7, its IRS 
spectrum is very similar to that of the composite spectrum of 
SMGs presented in Pope et al. (2008), and for which they report 
T9.7 ~ 1 in excellent agreement with what we expected based on 
this source's S 100/S 24 colour. We should also note that GN26 
was one of the 13 sources used for the construction of the com- 
posite spectrum. Taken together, the evidence suggests that the 
red S 100/ S 24 colour of our sample is caused by the existence of 
a moderate/strong silicate absorption feature at 9.7 /mi, that en- 
ters the 24 yum band at these redshifts. We therefore choose to 
characterize these sources as silicate-absorbed galaxies. 

4.3. The far-IR properties of the high-z sample 

In the pre-Herschel era the far-IR properties of large samples 
of galaxies were derived based on large extrapolations of their 
24 fim flux density using SED templates based on local IR- 
luminous galaxies. Recent studies using Herschel data have con- 
firmed the validity of these extrapolations, demonstrating that 
the 24yum flux density is a good proxy of the total L IR , at least 
up to z ~ 1.5 (Elbaz et al. 2010). On the other hand, here we 
have found a population of high-z sources (1.0 < z < 1.7), that 
exhibit atypically red S 100/S 24 colours. For these sources, we 
expect that L IR estimates based on their S24, and using average 
template SEDs, such as CE01 and DH02, would be severely un- 
derestimated. 

In Fig. 12 we plot the S m /S 24 colour of all GOODS-H 
galaxies as a function of the ratio between the Herschel based 
Lir and the one derived using only the 24//m flux density. 
Although for the majority of the sources the two Li R estimates 
are in good agreement (within a factor of ~ 1.5), for the high-z 
sources with S 100/ S 24 > 43 the 24 jim flux density would un- 
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Fig. 11. Full SED of GN26, a source at z = 1.21 with S m /S 2 4 ~ 70. The 
far-IR part of the SED is derived from Herschel data and is overlaid with 
the best fit CE01 model. The mid-IR part of the SED is the observed 
IRS spectrum from Pope et al. (2008), that is also presented separately 
in the inset figure. We also show the optical/near-IR part of the SED 
along with the best fit BC03 model. 
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indicate the area where the two Lj R estimates agree within a factor of 
1.5. 



derestimate the true L IR on average by a factor of > 3. This is 
not the case for the low-z sample where the two estimates are in 
better agreement. For the high-z sample there is a clear trend be- 
tween the S 100/ S 24 colour and the ratio of the two Lir estimates, 
in the sense that the true L IR is progressively underestimated for 
sources with redder S 100/^24 colours, depicting the limited va- 
riety of SEDs used for the derivation of the Lir. This indicates 
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Fig. 13. Left: The SFR - M* relation for sources with S 100/^24 > 43 at 1.1 < z < 1.7 (black circles). Orange dots indicate sources with spectroscopic 
redshift. Red cross corresponds to the ULIRG dropout that we presented in Section 3.3. The solid green line depicts the SFR-M* correlation at z 
~ 1.4 (Elbaz et al. 201 1), while the dashed lines indicate its dispersion. Right: Specific star formation rate (sSFR), as a function of redshift, for the 
same sample as in the left panel, as well as for the whole GOODS-H sample (grey dots). The gray thick line denotes the evolution of the sSFR with 
redshift, as derived by individual detections and stacking analysis of the GOODS-H sample by Elbaz et al. (2011). The blue shaded area along 
with the dashed gray lines which correspond to the dispersion of the evolution (~ 0.3 dex), indicate the region of main sequence galaxies. Sources 
with S ioo/5 24 > 43 tend to be found above the blue shaded area, populating the starburst region. The green text indicates the position of Arp 220 
and M82 (shifted to z = 0. 1 for clarity) 



that although the average template SEDs are representative for 
the bulk of the galaxy population at z < 1.5, this is not the case 
for a population at 1.0 < z < 1.7 with very red S 100/^24 colours 
(> 43 for this study). 

4.4. The starburst nature of the high-z sample 

As the silicate absorption feature merely requires a mass of 
warm dust obscured by a significant column of cooler dust, it 
does not provide any insight into the mechanism that is heat- 
ing the warm dust. Hence, it is difficult to establish a correlation 
between the strength of this feature and the source that powers 
the mid-IR emission of infrared luminous galaxies, as it could 
equally be produced by a deeply buried AGN or a compact star- 
burst (Farrah et al. 2008, Imanishi et al. 2009, Armus et al. 2007). 

We have already reported than none of our sources show di- 
rect signs of AGN activity, since none of them are either (a) de- 
tected in X-rays in the 2 Msec Chandra observations (Alexander 
et al. 2003), (b) satisfies the criteria of power-law AGN or (c) 
has high excitation lines in their optical spectra (where avail- 
able). A deeply obscured AGN though, cannot be ruled out. 
To investigate this, we stack the Chandra X-ray data on these 
sources. We find a strong detection (7<x) in the 0.5-2 keV band 
corresponding to L2-10 = 3 x 10 41 erg sec" 1 and no detection in 
the 2-7keVband (3cr upper limit 1.28 xlO" 17 erg cm -2 s" 1 ). The 
derived upper limit in the hard band along the stacking results in 
the soft band suggest that the sources are dominated by star for- 
mation (Nandra et al. 2002, Lehmer et al. 2008). Furthermore, in 
most cases the radio based L IR is in good agreement with that de- 
rived by Herschel, again indicating that those sources are domi- 
nated by a nuclear starburst, rather than an AGN. 

It has recently been shown that normal star-forming galax- 
ies exhibit a correlation between their SFR and stellar mass at 



any given redshift. This correlation was first found among z ~0 
galaxies (Brinchmann et al. 2004) and was subsequently ex- 
tended to higher redshifts, z ~ 1 Elbaz et al. (2007), Noeske 
et al. (2007), z ~ 2 Daddi et al. (2007), Pannella et al. (2009), 
z ~ 3 Magdis et al. (2010a,b,c) and z ~ 4 Daddi et al. (2009). 
It has also been shown that star forming galaxies that do not 
follow this correlation tend to undergo a rapid starburst phase 
and have more compact geometries (Elbaz et al. 2011). For in- 
stance, at z ~ and z ~ 2, respectively, local ULIRGs and 
SMGs have SFRs that greatly exceed the SFR-M* correlation. 
Both classes of objects are thought to host compact starbursts 
(Daddi et al. 2010, Tacconi et al. 2010). It would therefore be 
interesting to investigate the position of the sources with red 
S 100/ S 24 colour in the SFR - M* diagram. Using the Kennicutt 
1998 relation and a Salpeter IMF, we convert the Herschel based 
L IR to SFR for the sources with S 100/ S 24 > 43 and 1 . 1 < z < 
1.7 and we plot the derived SFR versus the stellar mass of the 
galaxies (Fig. 13, left). We also overplot the SFR-M* correla- 
tion at the median redshift of the sample i.e., z = 1.4 (Elbaz et 
al. 2011). Clearly the sources are off the correlation, exhibiting 
enhanced star formation activity for their stellar mass. Similarly 
the ULIRG dropout source that we presented in section 3.3 is 
also off the correlation. To eliminate the effect of the evolution 
of the SFR-M* with redshift, we also show the specific star 
formation rate (sSFR), defined as SFR/M*, as function of red- 
shift for sources with S100/S24 > 43, as well as for the whole 
GOODS-H sample (Fig. 13, right). We also overplot the evolu- 
tion of sSFR with time as derived by Elbaz et al. (2011), based 
both on detected sources and stacking analysis of the GOODS- 
H sample. According to Elbaz et al. (2011), sources within the 
blue shaded area are main sequence, normal star forming galax- 
ies galaxies while sources above it, are considered to undergo a 
starburst phase. We find that sources with S 100/^24 > 43 tend to 
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Fig. 14. The fraction of sources expected to be missed in the 24 /urn band at 0.2 < z < 0.6 (left) and 1.0 < z < 1.7 (right) as a function of the 24 /urn 
vs 100 yL/m (top) and the 24 yi/m vs 160 yi/m (bottom) depths. Black circles indicate the depths of PACS and MIPS24 yi/m surveys for some of the most 
important cosmological fields, i.e EGS, Lockman and COSMOS. We note that the region where the constant fraction isocontours run horizontally, 
cannot be constraint by the data. The adopted detection limits (5 <x) at 24 yi/m are 80 yt/Jy for COSMOS and 50 yt/Jy for EGS and Lockman. Similarly 
the 3 cr detection limits at 100 ^m are 5.0-, 3.8- and 3.6mJy (COSMOS, EGS and Lockman), while at 160 yt/m are 10.2-, 8.6- and 7.5mJy. 



have higher specific star formation rates when compared to main 
sequence star forming galaxies at this redshift range, and popu- 
late the starbursts region. This result is coherent with the picture 
where these sources are compact starbursts with high sSFR and 
where (like local ULIRGs) a strong silicate absorption feature at 
9.7 fim is present in their mid-IR spectra. 



5. Discussion 

Previous studies predicted that the number of such silicate break 
sources could reach ~ 900 to 1500 sources per square degree, 
depending on the assumed model (e.g. Tagagi & Pearson 2005). 
Furthermore, Kasliwal et al. (2005), based on 16yum IRS data, 
reported that such sources account for more than half of all 



galaxies at z ~ 1-2 predicted by various models. The samples of 
silicate-break galaxies and MIPS dropout sources that we have 
found here can place strong limits on the number of 1.0 < z < 
2.0 infrared luminous galaxies with similar properties. 

In the 250 arcmin 2 covered by PACS in the two GOODS 
fields, we have detected -30 (7) silicate absorbed candidates at 
1.0 < z < 2.0 with S m /S 2 4 > 43 and L lR > 10 11 L Q (L IR > 
10 12 L©). Assuming that all sources with S 100/S 24 > 43 at this 
redshift are silicate absorbed starbursts, and neglecting the ef- 
fects of cosmic variance, this implies a surface density of ~ 540 
(220) sources deg" 2 , significantly lower than previous predic- 
tions. Using the co-moving volume within 1 < z < 2, this im- 
plies a space density of Si absorbed ULIRGs in this redshift bin 
of ~ 2.0 (± 0.3) x 10" 5 Mpc" 3 . Comparing the numbers with 
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Fig. 15. The fraction of sources expected to be missed by the 24 iim band at all redshifts as a function of the 24//m vs 250 iim and the 24 jjm vs 
to 350 yt/m depth. Black circles indicate the depths of PACS and MIPS24 yi/m surveys for some of the most important cosmological fields, i.e EGS, 
SWIRE and COSMOS. 



the rest high-z sample of GOODS-H, we find that these sources 
account for the ~ 8% and 16% of the Lir > 10 11 L Q and Lir > 
10 12 L sources at 1.0 < z < 2.0. We note though, that our sam- 
ple is not complete. As illustrated in Fig. 9, sources with sili- 
cate absorption feature but cold T^ ust have S 100/^24 < 43 and are 
missed from the selection. We therefore conclude that this esti- 
mate is a rather conservative lower limit. Moving to the dropout 
sample, we have identified 11 sources with Lir > 10 11 L©, ac- 
counting for approximately 1-2% of the population of 1 < z < 
2 infrared luminous galaxies in GOODS-H as a whole. We note 
that by comparing to the whole GOODS-H sample in the same 
redshift range, all effects of incompleteness are taken into ac- 
count. Finally, these number should be treated as upper limits, 
given that some of these could be spurious detections, as dis- 
cussed in Section 2. 

These estimates confirm that for the GOODS surveys, 
24 [im observations recover the vast majority of z < 2 sources 
and do not introduce a strong selection bias. Furthermore, we 
confirm that using the 24 jim catalogues to define priors for the 
extraction of PACS sources should only miss a small fraction 
of high-z sources. Of course the ratio between the number of 
sources missed by the 24//m band and the total number of 
sources in the PACS bands strongly depends on the depth of the 
MIPS 24 yum and PACS observations. This is illustrated in Fig. 
14, where we show the fraction of the expected dropouts at 0.2 
< z < 0.6 and 1.0 < z < 1.7 as a function of the 24- to 100 yum 
and the 24- to 160/mi depth. These estimates are solely based 
on observations (i.e. on the GOODS-H sample), without any 
further assumptions. We simply calculate the ratio of /drop- 



number of sources with S 



100 



> a & S 2 4 < P 



number of sources with S100 > ol & S24 > J3 

where a and ft are free parameters. These diagrams can be 
used as diagnostic of the fraction of sources that will be missed 
by a source extraction method based on 24 fim priors for sev- 
eral extragalactic Herschel surveys, e.g. PEP (Lutz et al. 2011), 
HerMES (Oliver et al. in prep). We conclude that for the ma- 



jor extragalactic surveys, there is not a large population of sil- 
icate break galaxies that would have been undetected in 24 /mi 
Spitzer data as some authors have previously suggested. This dis- 
crepancy could be indicative of an evolution of the strength of 
the silicate absorption feature with time. If the silicate features 
were weaker in early galaxies due to more extended star forming 
regions, then theoretical models based on local templates would 
grossly overestimate the number density of such sources. This 
scenario is in line with recent findings, where the star formation 
activity of the majority of high-z galaxies (including those with 
high Lir), is similar to that of local spiral galaxies (e.g. Daddi et 
al. 2010, Genzeletal. 2010). 

Finally, although in this work we have mainly focused on 
PACS data, it is worth attempting to quantify the fraction of 
sources that would be missed by a 24 jim prior based source ex- 
traction in the SPIRE bands. The philosophy behind Fig. 15 is 
identical to that of Fig. 14, but this time considering the whole 
redshift range of the sources and the S250 (left) an d S350 (right) 
flux density limits. We note that these detection limits, indicate 
the expected instrumental noise and do not take into account the 
confusion noise, which is dominant in the SPIRE bands. 



6. Conclusions 

We have presented a study dedicated to sources that exhibit atyp- 
ical S ioo/S 24 colours using the deepest PACS data to date, ob- 
tained as part of the GOODS-H program. By performing blind 
source extraction we searched for sources that are bright in the 
far-IR but undetected at 24 /mi, i.e. for 24 fim dropout galaxies. 
Then we investigated the properties of sources in the GOODS-H 
sample with red S 100/S 24 colours searching for the population of 
silicate-break galaxies. The main results are summarized below: 

- We have identified 21 PACS sources that are undetected at 
24yum (down to a 3<x detection limit of ~ 20 yuJy). These 
24 /mi dropout sources are found to have a bimodal redshift 
distribution, with peaks centred at z ~ 0.4 and ~ 1.3, and 
are expected to exhibit strong silicate absorption features, 
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responsible for their depressed 24 jim emission. Among the 
sources in higher redshift we identify 10 LIRGs and one 
ULIRG at z =1.68. This enables us to place upper limits in 
the fraction of LIRGs/ULIRGs that are missed by 24 jim sur- 
veys. 

- The vast majority of Herschel PACS sources are detected at 
24 yum, indicating that a prior-based source extraction based 
on the 24 yum emission of the galaxies suffers only very mod- 
est incompleteness, with MIPS dropout sources accounting 
only for ~ 2% of the infrared luminous population in the 
GOODS fields. Although this fraction is negligible for the 
GOODS surveys, 24 jim dropouts may be a concern for other 
Herschel extragalactic surveys with shallower 24 jim data. 

- Based on the mid- and far-IR colours of sources in the 
GOODS-H sample, we demonstrated that sources with 
S100/S24 > 43 and Sig/S$ > 4 are located in a narrow red- 
shift bin, 1.0 < z < 1.7. Furthermore, we provided evidence 
that sources selected in this manner are starburst- dominated 
and with compact geometries. Similarly to the dropouts, the 
red S ioo/S 24 colours of these sources are attributed to the 
9.7 fim silicate absorption feature in their mid-IR spectra that 
enters into the 24 jim band. We characterize them as silicate- 
absorbed galaxies. 

- The infrared luminosity of these silicate-absorbed galaxies, 
when derived based on their monochromatic 24 jim flux den- 
sity, is on average underestimated by a factor ~ 3. They ac- 
count for about 16% of the ULIRGs in the GOODS fields 
with a space density of 2.0 x 10" 5 Mpc -3 . 

- We provide diagnostic diagrams to estimate the fraction of 
sources expected to be missed in the 24 yum band for several 
Herschel extragalactic surveys, and predict that for most of 
them that fraction is less that 10%. 
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